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Abstract. We have developed a new ISOCAM data re- 
duction technique based on wavelet analysis, especially 
designed for the detection of faint sources in mid-infrared 
surveys. This method, the Pattern REcognition Technique 
for Isocam data (PRETI) has been used to reduce the ob- 
servations of the Hubble Deep Field (HDF) and flanking 
fields with ISOCAM at 6.75 (LW2) and 15 /jm (LW3) 
Rowan-Robinson et al. (|1997|). Simulations of ISOCAM 



greater uncertainties, given the small number of detec- 
tions. 



data allow us to test the photometric accuracy and com- 
pleteness of the reduction. According to these simulations, 
the PRETI source list is 95% complete in the 15 /im band 
at 200 /iJy and in the 6.75 /im band at 65 /xJy, using 
detection thresholds which minimize the number of false 
detections. We detect 49 objects in the ISO-HDF at high 
confidence secure level, 42 in the LW3 filter, 3 in the LW2 
filter, and 4 in both filters. An additional, less secure, list 
of 100 sources is presented, of which 89 are detected at 
15 /im only, 7 at 6.75 only and 4 in both filters. All 
ISO-HDF objects detected in the HDF itself have optical 
or infrared counterparts, except for one from the addi- 
tional list. All except one of the r adio sources detected in 
the field by Fomalont et al. ( 1997 ) are detected with ISO- 
CAM. Using a precise correction for the field of view dis- 
tortion of ISOCAM allows us to separate blended sources. 
This, together with the fact that PRETI allows to correct 
data on the tail of cosmic rays glitches, lead us to pro- 
duce deeper source lists than previous authors. Our list of 
bright sources agree with those of Desert et al. (1998) in 
both filters, and with those of Goldschmidt et al. (1997) in 
the LW3 filter, with systematic difference in photometry. 
Number counts derived from our results show an excess by 
a factor of 10 with respect to the prediction of a no evo- 
lution model (Franceschini 1998) in the LW3 band. On 



the contrary, the number of sources in the LW2 band is 
compatible with the prediction of such a model, but with 
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1. Introduction 



The observation with the Hubble Space Telesco pe (H ST) 
of the Hubble Deep Field (HDF) (Williams et al. |l996|) has 
yielded numerous new results in observational cosmology. 
This is due to its depth and spatial resolution in the near- 
UV and optical, allowing the computation of photometric 
redshifts of the observed galaxies as well as the determina- 
tion of their morphological type. With the ground -based 
follow-up observations, in the near IR (Hogg et al. 1997a, 
Connolly et al 1997) for example, or at radio wavelength 



(Fomalont et al. 1997), as well as spectroscopic measure- 
ments (Steidel et al. |l996j , Cohen et al. |l996j ), this field 



now constitutes a very important database for the study 
of galaxy formation and evolution. Among the striking re- 
sults that have been already produced (see Ellis ( 1997| ) 
for a review), two have strong implications on the galaxy 
formation picture. First, the star formation rate (SFR) 
at high redshifts has been derived (Madau et al. 1996[ ), 
and together with that computed from Canada France 
Redshifts Survey (CFRS) data at lower redshifts (Lilly et 



al. 1996), it seems to indicate that the bulk of star for- 
mation arose below z ~ 2. Second, an increase of faint 
irregular galaxy count s around z « 1 has been reported 



(Abraham et al. 1996 ), as well as the lack of detection of 
quiescent evolving early types galaxies at high redshifts. 
These results tend to favor the hierarchical galaxy forma- 
tion picture. 

However, the interpretation of these results is now be- 
ing widely discussed. For instance, it has already been 
pointed out that the fainter objects seen in the HDF could 
merely be star forming regions of galaxies and not entire 
galaxies, because the distance between detected objects 
is small with respect to the size of a galaxy (Colley et 
al. 1996 ). This idea is supported by new simulations that 
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show that diffuse emission from evolved stellar popula- 
tions is below the detection limit of the HDF, while star 
forming regions with high UV surface brightness can be 



developed a new technique for ISOCAM data reduction, 
based on wavelet analysis (Starck et al. 1998) : the Pat- 



detected (Hibbard & Vacca 1997). For similar reasons, qui- 
escent evolving early-type galaxies at redshift of z ~ 3 are 



also undetectable (Maoz 1996). 

Most of the data that were used to derive star forma- 
tion rates are UV and optical data from the CFRS and 
the HST itself, which are sensitive to extinction. Indeed, it 
has been shown that star formation at high redshift is ob- 
scured, typically by ~ 2 to 3 mag at A = 1620A(Meurer et 
al. 1997). The SFR inferred might therefore change if the 



reprocessing of UV by dust into IR is taken into account. 
Knowledge of the dust content of galaxies is therefore rel- 
evant to assess the results obtained so far, both at z < 1 
with the CFRS and at higher redshifts. Additional infor- 
mation on this problem can be obtained by observing the 
HDF field in the infrared. 

Many models have been built to compute the in- 
frared spectra of various types of galaxies (AGN, star- 
bursts, early types, etc) and make predictions on the num- 
ber counts for surveys in the infrared (see, for example, 
Franceschini et al. ( |1991 ). Moreover, various groups have 
derived new IRAS 60 fim counts, going to the extreme lim- 
its of the sensitivity of this satellite. Their results show an 
excess of sources with respect to the predictions of a sce- 
nario with only passive luminosity evolution of galaxies 
(Bcrtin et al. 1997, Gregorich et al. 1995). However, the 



exact amount of evolution needed is still a matter of de- 
bate, because the number counts derived by the various 
teams are slightly different. This issue needs to be revisited 
on deeper samples. From this point of view, the HDF field 
is a great opportunity to test these models, first because 
their number count predictions can be checked against ob- 
servations, and second because all the information at other 
wavelengths is of great help to test the hypotheses of the 
various models. 

Since the HDF field is rather small (~ 5 arcmin 2 ) and 
very deep, infrared observations require a good spatial res- 
olution as well as a good sensitivity. This can be achieved 
by the mid-IR imager ISOCAM of the Infrared Space Ob- 



servatory (ISO) (Kessler et al. 1996), which has proven to 



be 1000 times more sensitive than IRA S, wit h a 50 times 
finer spatial resolution (Cesarsky et al. 1996 ). Part of the 



Director's Time was dedicated to complete these obser- 
vations, at 6.75/xm and 15/xm, for a total amount of 12.5 
hours, with M. Rowan-Robinson as Principal Investiga- 
tor (Serjeant et al . citeSerjeant, Goldschmidt et al. 1997 , 
Oliver et al. |l997i Mann et al. |l997|, Rowan-Robinson et 



al. 1997. The data of the ISO-HDF are now in the public 



domain, and we present here the results of our processing. 

The ISOCAM observation of the HDF, hereafter ISO- 
HDF, is at the moment the deepest survey obtained in 
the 15 /xm band. It is also deep in the 6.75 fxm band, al- 
though a deeper integratio n on t he Lockman Hole was ob- 
tained by Taniguchi et al. ( 1997 ) with this filter. We have 



tern REcognition Technique for ISOCAM data (PRETI). 
Compared to data that we have at hand from less deep 
surveys in the Lockman Hole, from the Guaranteed Time 
program by Cesarsky et al. (1996), and in the shallow 



survey of the ELAIS program, this set has a much higher 
redundancy, and thus is useful for finding the ultimate ca- 
pacities of our data reduction method. In addition, since 
the ISO-HDF has been analyz ed by other teams s uch a s 
Rowan- Robinson et al. ( 1997 ) and Desert et al. ( 1998 ), 
this observation allows us to compare the results from our 
method to those of more classical ones. All these methods 
aim at extracting very faint sources in ISOCAM data, and 
independent reduction techniques improve the quality of 
the produced catalogs. In section 2, the characteristics of 
the emission expected from galaxies in the MIR are out- 
lined, and compared to the bandpasses of the ISOCAM 
filters that were used for the HDF observation. In section 
3 we present the data acquisition and stress the critical 
steps of the reduction process. The way in which these 
difficulties are handled by the PRETI method is described 
in section 4. The efficiency of the method for source detec- 
tion and its photometric accuracy was also tested on sim- 
ulated data, which are presented in a section 5. Finally, 
we give our source catalog and number counts, compare 
them with those of other groups, and discuss the results. 
A more detailed interpretation of the results, source by 
source, will be presented in a forthcoming paper (Elbaz et 
al., in prep.). 

2. MIR emission in ISOCAM filters 

The ISO observations were performed in two broad band 
filters, LW2 (from 5.0 (im to 8.5 /im) and LW3 (from 12.0 
//m to 18.0 /im). Figures and ^ present synthetic spectra 
of the mid IR emission of galactic components between 1 
and 20 /im. Figure [l] is a spectrum of a starburst galaxy, 
taken from Moorwood et al. (1996) for the 1 to 5 /im part 
and from Tran( |1998[ ) for the 5 to 18 fim. part. Unidenti- 
fied infrared bands (UIB) and ionization lines of Neon are 
clearly visible, as well as a hot dust continuum and silicate 
absorption. In this case, the derived temperature for the 
hot dust is about 200 K, using the model from Desert et 
al. (1990), and the extinction due to dust is about 3 mag, 



representative of a starburst region in a galaxy. The dot- 
ted line is a spectrum dominated by very hot dust heated 
at 1000 K, that could be typical of the mid-IR emission 
in t he cen tral core of an AGN, as observed by Aussel et 
al. ( |1998| ). Figure | depicts emission dominated by old 
population stars as modeled by Bressan et al. ( |1998 ), and 
is thus relevant for early- types galaxies. 

The mid-IR (hereafter MIR) colours of galaxies are 
not simple, because the observed emission of a galaxy is 
a mixture of all the components shown on Fig. |l| and Fig. 
||. According to the type of the object (early-type, AGN, 



H. Aussel et al.: ISOCAM observations of the Hubble Deep Field reduced with the PRETI method. 



3 



1 .00 



UW2 Z =^ T LW3 2=1 T 



^LW2 z = 0.5^ j LW3 z = 0.5 j 



LW2 z = 



LW3 z = 



11.3 



12.7/NeII 



NeHI 



0.10 — 



0.01 




AGN' 



HOT DUST (VSG) 
T ~ 200 K 



1 

Wavelength ((in) 



20 



Fig. 1. Sketch of the mid IR emission of a starburst galaxy (after an original sketch of A. Jones) and of an AGN, 
compared with ISOCAM LW2 and LW3 filters bandpasses at various redshifts. 

Solid line: total MIR emission of a starburst galaxy. The 5/im - 18/im part is a an ISOCAM CVF observation showing 
all the typical UIB bands and ionization lines (spectra kindly provided by D. Tran). The 2/xm - 5/im part was taken 
from a SWS spectrum of the Circinus galaxy (Moorwood et al. 1996 ), and scaled to the CVF spectrum. 
Dashed line: contri bution to the MIR emission of the hot dust. This dust is constituted of very small grains (VSG) 
(Desert et al. 1990|) heated to a few hundred degrees. 

Dotted-dashed line: Possible contribution of a continuum associated to the carriers of UIB. 

Dotted line: Hot dust (1000 K) MIR emission of an AGN central part, as observed by Aussel et al. ( |1998| ). 



starburst or normal galaxy), it is expected that a given 
type of emission will dominate (e.g. stellar, very hot dust 
or UIB). T his is not always the case, as shown by Mad- 
den et al. ( 199§| ) on a sample of elliptical galaxies : most 
exhibit stellar colours in the MIR but some present MIR 
spectral energy distributions (SED) typical of UIB emis- 
sion. Thus, the relative fluxes in figs. |l| and || are not to 
be compared directly. 

In the rest frame, the LW2 filter is dominated by UIB 
features while the emission in the LW3 filter is due to 
the hot dust continuum, ionization lines of neon and the 
12.7 /jm UIB plus the associated continuum. This picture 
changes dramatically when the redshift rises toward z=l : 



the contribution of the stellar continuum, especially from 
old population stars, overtakes UIB features in the LW2 
band, which in turn are shifted to the LW3 band, together 
with the silicate absorption band at 9.7 /zm. The wave- 
length coverage for both filters is shown on both Fig. [I] 
and Fig. ||for redshifts of 0, 0.5, and 1. This clearly indi- 
cates that the mid-IR observations can only be interpreted 
with the help of data at other wavelengths to disentangle 
the various contributions. This is possible on the HDF 
field, because UV, optical and near IR, as well as radio 
observations are available. 
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Fig. 2. Sketch of the mid IR emission of an old stellar population, compared with ISOCAM LW2 and LW3 filters 
bandpasses at various redshifts. 

Solid line : model of the MIR emission o f a 5 Gyr old stellar simple population, where dust envelopes of Mira and 
OH/IR stars are neglected (Bressan et al. 199§| ). 

Dashed line : model of the MIR emission of a 5 Gyr old stellar simple populati on, for which the dust envelopes have 
been taken into account for the AGB phase, using silicate grains (Bressan et al. 199S). 

The arbitrary units used for fluxes cannot be compared with those of Fig. |l|, the relative importance of these emissions 
in various galaxies varies by factor of order 100. 



3. Data acquisition 

The HDF field is well suited for mid-infrared deep observa- 
tions. It was selected for its low optical extinction, low HI 
column density and, most important fact for ISOCAM ob- 
servations, low c irrus emission in the IRAS 100 /jm maps 
(Williams et al. 1996 ). Moreover, it sits at high galactic 
latitude, and contains therefore only a few stars. Simi- 
lar criteria were used to select fields for other ISOCAM 
deep s urvey s such as the IDSPCO program (Cesarsky & 
El baz |1996|) or the DEEPPGPQ program (Taniguchi et 
al. 



1997|) . 



Table [l] describes the observation parameters for both 
filters. Each consist of three 8x8 rasters, which are partly 
overlapping. In both cases, the area covered is wider that 



the HDF field, and even contains a large part of the flank- 
ing field for the LW3 observation, which was done with a 
larger pixel field of view (PFOV) and a greater raster step 
than the LW2 observation. If the detector had been ideal, 
i. e. with no memory effect and insensitive to cosmic rays, 
a very good sensitivity would be reached : 



, a . , 

pixels pixel 
y exp 

where a P i xe i is the noise per readout for each ISOCAM 
pixel (it is therefore equal to the readout noise plus the 
photon noise, which is dominated in the case of faint 
sources by the background emission) , and N exp is the 
number of exposures. One can easily compute the ideal 
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Table 1. ISO-HDF observation parameters. The main pa- 
rameters for each raster and for the whole observation are 
listed 



For the ISO-HDF data, as for all other data, the process 
goes through the following steps : 



Parameter 


LW2 (6.75 Aim) 


LW3 (15 Mm) 


Raster 


Integration time 


10 s 


5 s 


Detector Gain 


2 


2 


Number of exposures 


9 or 10 


19 or 20 


Pixel field of view 


3 " 


6 " 


Number of pointings 


8x8 


8x8 


Step between two 


5" 


9 " 


pointings 




Observation 


Number of rasters 


3 


3 


Total area covered 


10.45' 2 


24.3' 2 


Area of maximum 


0.98' 2 




sensitivity 


Sensitivity (Jy/" 2 ) 


a < 0.2 


a < 0.34 


Exposure time 


> 1.94 hrs 


> 2.25 hrs 


Area of mean sensi- 
tivity 


7.25' 2 


21. 5' 2 


Sensitivity (Jy/" 2 ) 


a < 0.5 


a < 1.37 


Exposure time 


> 0.31 hrs 


> 0.17 hrs 



1. 



translate the raw data (CIER and IIPH files) in a 
raster structure containing the images, the ISOCAM 
configuration and the pointing information, using Cam 
Interactive Analysis (CIA). 

convert the images in normalized units of ADU/G/s 
and subtract the dark current, with CIA. 
correct the data for short time cosmic rays and mask 
the affected readouts. 

clean the data for long term effect of cosmic rays. 

subtract the baseline of each pixel. 

compute the reduced images per raster pointings. 

derive a flat-field from the baseline and correct the 

reduced images accordingly. 

project each raster pointing on the sky. 

perform a spatial detection on the final image. 



The general method is described in Starck et al. (1998) 
and does not need to be detailed here. However, we discuss 
here some of the points above, either because they are very 
important, or because they need to be specifically adapted 
to each observation. 



4-.1. Long term effect of cosmic rays. 



sensitivity of the ISO-HDF if one computes a P i xe i for each 
pixel : this is easily done by fitting a Gaussian distribution 
on the histogram of the each pixel, after having filtered 
it of all the low frequency components, because all the 
non ideal behaviour of ISOCAM occurs mostly at low fre- 
quency. At a 5cr level, one could reach 15 [±iy in the LW3 
band and 5 /iJy in the LW2 band. The observation was 
designed to be as deep as possible, and with a good spatial 
resolution in order to be able to identify with confidence 
the optical counterparts of ISOCAM sources. The observa- 
tions were therefore made using the microscan mode, that 
is by moving the camera by fraction of pixels between two 
raster steps. With a step of 5 arcsec and a PFOV of 3 arc- 
sec, the final resolution of the ISOCAM map should reach 
1 arcsec in LW2, and 3 arcsec in LW3. 

Unfortunately, the long wavelength detector (LW) of 
ISOCAM is very sensitive to cosmic rays due to its thick- 
ness, and shows a transient behaviour. These two main 
drawbacks alter its efficiency in detecting faint sources. 
This is the reason why the PRETI method was developed. 



4. Data reduction with the PRETI method 

The two main difficulties for the detection of faint sources 
with the ISOCAM LW channel are mostly the cosmic ray 

19981 ) . 



effects and the transient behaviour (Starck et al. 
The PRETI method has been designed to overcome these 
two problems. The main idea of this method is to clean the 
cube of ISOCAM data as well as possible, before recon- 
structing the images and performing the source detection. 



As explained in Starck et al. (1998), the impact of cos- 
mic rays, apart from the high frequency response they 
induce in the LW detector, are sometimes followed by a 
slow decrease of the signal ('faders') or by a depletion in 
the detector gain ('dippers'). This behaviour is attributed 
to proton or a particle impacts on the detector, while sim- 
ple 'glitches' are due to effects of the cosmic ray electrons. 
One has to be very careful in cleaning them, because they 
will produce false detection on the final map if they are 
not well removed. 

This is especially true for the ISO-HDF observations, 
because only a small number of readouts are done per sky 
position, and because of the microscan observation mode. 
In this observation, one third of the 'faders' and 'dippers' 
lasts more than 60 readouts. This means that they last for 
more than 6 pointings (LW2) or 3 pointings (LW3). But 
each cosmic ray hits on average 4.5 adjacent pixels, that 
will more or less have the same behaviour. In a microscan, 
each position on the sky is observed by two adjacent pixels 
during two successive pointings. Figure ^ illustrates such 
an event : if the 'fader' effect is not corrected or masked 
out, it produces a spurious source when averaging all the 
readouts of a same pointing, because the flux appears to be 
above the 5 a noise level obtained on individual pointings. 
This example is taken from the LW3 observation : the 
two first pointings affected by the cosmic rays are less 
distant that the PSF FWHM. Since adjacent pixels of the 
one in the figure present the same behaviour, a spurious 
source will be 'detected', because faders will be co-added 
together. 
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Fig. 3. Effect of a cosmic ray impact on the ISOCAM LW detector. The 'glitch' occurred on readout 9, and produced 
a 'fader' effect. The solid horizontal line at 28.85 ADU is the average background level. The solid horizontal line at 
29.92 is the 5 a level of the noise. The dash-3 dots horizontal line is the average 5 a noise level after averaging all the 
readouts of a pointing. The dashed-dotted horizontal line is the 5 a level on the final raster map. The limits of each 
ISOCAM pointings are represented by dotted vertical lines. This 'fader' last for 3 pointings. 

solid line : pixel response versus readout number, after correction for the short term effect of cosmic rays. Frames 9 to 
15 and 26 to 30 have been masked out. 

dashed line : final pixel response after correction for the long term effect of cosmic rays. 



This effect is even more important in very deep obser- 
vations where many raster positions are coadded together : 
'faders' and 'dippers' are detected in the time history of 
the pixel, thus against a rather high noise level (see Fig. ||). 
On the contrary, sources will be detected on the co-added 
map, with a very low noise level that has been represented 
on the same figure. It may happen that the tail of cosmic 
ray effects, that were undetected in the temporal history, 
are coadded to themselves and therefore become visible 
on the final image as spurious structure. 

In order to confront these problems for the HDF sur- 
vey, we put very strong constraints on the deglitching pro- 
cess, that is : 



1. The detection of structures in the temporal history of 
the pixel is done against 4 a of readout and photon 
noise, a lower level than for the usual PRETI method. 

2. Each positive structure whose maximum is expressed 
at the same scale as a raster pointing is treated as a 
cosmic ray if : 

a) a cosmic ray has been detected at its beginning 
(i.e. the two first readouts have been masked by the 
short deglitching process). 

b) more than 50% of the readouts are already 
masked by the short deglitching process 

This intends to avoid as much as possible false detec- 
tions. There is the possibility that we might clean out a 
source, but this can be checked by the following method : 
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for each data set that we treat, we produce also raster 
maps where all the long term effects of cosmic rays are 
masked out instead of being corrected. Of course, this 
leads to far less sensitive maps since 20 to 40 % of the 
readouts are then eliminated. But this allows to check that 
the photometry of the remaining sources is consistent be- 
tween the two reductions : this ensures that no real source 
was 'cleaned out' during the reduction. Indeed, in the end, 
we never find discrepencies in the flux of a given source 
by factors greater than 15%. 



4-2. Baseline subtraction. 

In order to correct for the long transient behaviour that 
takes place at the beginning of each observation (due to 
the change of the ISOCAM configuration), and to clean 
the data of the possible drifts of the detector, the PRETI 
method suppresses the baseline of each pixel's time his- 
tory. This is computed at larger temporal scales of the 
signal, at least larger than the scale of one raster pointing 
(Starck et al. 



1998) 



Because of the microscan observation, one has to ex- 
tract the baseline at least two scales further for the ISO- 
HDF observation. Since the step between two raster point- 
ings is about the same size as the FWHM of the ISOCAM 
PSF (9 arcsec at 15 /xm and 4.5 at 6.75 /im), each source is 
seen by the same pixel during more than a raster pointing. 
In order to avoid to take some flux from the source into 
the baseline, it is computed two scales further as usual. 

It is worth noting that the subtraction of the baseline 
frees the result of the reduction from the possible errors 
in the dark subtraction. Final images that are produced 
by the PRETI method are always on a zero background. 
Thus, such a method is only usable for observation on 
constant background, with no structure. This is the case 
for the zodiacal light emission, on an area as small as the 
HDF. 

Nevertheless, it is possible to construct a background 
map of the observation. This map can be used to check 
the photometry by comparing the observed background 
with predicted values of models of zodiacal light. 



4-3. Projection onto the sky. 

Once each raster position has been averaged and flat- 
fielded, it is projected onto a sky map (RA, DEC, epoch 
2000), whose resolution is computed according to the step 
size between raster pointings. The algorithm being used is 
derived from the classic 'shift-and-add' method. The in- 
tersection area of each CAM pixel with the sky map is 
computed, and added to it with a weight that is equal to 
this area times the square root of the number of pixels 
that were coadded to give the raster pointing image. That 




Fig. 4. ISOCAM distortion pattern for the LW10 6" 
PFOV configuration. Each vector ends where the center 
of a pixel should fall if they were no distortion and starts 
from its observed position. The length of the vectors are 
at the scale of the plot : the effect is sometimes greater 
than the size of one pixel. 



E 



pointings ' 



E 



pointings 



Assuming Gaussian distributions for pixels, one has for 
the noise map : 



E 



pointings ^{x i y 1 i i j)'^'^3^ 



h3 



S 2 N- ■ 

pointings (x,y,i,j) % >3 



where R x . v is the value of the final raster map at (x,y), 
S( x ,y,i,j) is the intercepted surface between ISOCAM pixel 
(i,j) and raster map pixel (x,y) and Nij is the number of 
readouts that where averaged together to produce Iij, 
the image of the raster pointing. We choose to use y/Nij 
as weight rather than the R.M.S. computed during the 
coaddition because this is very poorly estimated (less than 
10 readouts for the LW2 observation for example). Since 
the ideal detector should follow a Gaussian behaviour, its 
noise should vary as yJN^j , which is a more reliable esti- 
mator. It is straightforward from the computation of a x . y 
that the noise map against which we do our detection is 
highly inhomogeneous. 

The projection of each raster pointing takes into ac- 
count the field distortion of ISOCAM. It has been mea- 



sured by our team (Aussel [1999D , using the ISOCAM PSF 
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Fig. 5. Zoom on a part of the ISOCAM Hubble Deep Field observation at 15/im. The same processing was applied to 
the two images, except for the correction of the distortion. Cuts on images are at the same level. Pixel size is 2" . 



model from Okumura ( 1997 ). Followi ng the work done on 
the HST WFPC by Holtzman et al. ( |l995| ), we fit on our 
measurement a general polynomial of degree 3, that allows 
to compute for any pixel on the detector its position on 
the sky : 



ao + a\x + a-iy + a^x + a±xy - 
clqX 3 + a-jx 2 y + a%xy 2 + agy 



a 5 y 



y c = b + b\x + b 2 y + by,x 2 + b 4 xy + b^y 2 + 
b 6 x 3 + b 7 x 2 y + b s xy 2 + b 9 y 3 

Figure [| shows a map of the distortion of the LW chan- 
nel of ISOCAM with the 6" PFOV and the LW3 filter, 
where each vector ends where the center of a pixel should 
fall if they were no distortion, and starts from its actual 
position. The length of the vectors are at the scale of the 
plot : at the lower corners of the array (lines to 5), 
the effect is bigger than one pixel. If this effect were not 
taken into account, the sources would be smeared when 
coadding, and the final astrometry would be poor, reach- 
ing up to a 6" error. This is illustrated by Figure ||, where 
a zoom on the outskirts of one ISO-HDF raster is dis- 
played without and with correction of distortion. On the 
borders and corners of the raster, the gain of signal over 
noise ratio is about 50 %. 

The measurements at various wavelengths have shown 
that the distortion is a chromatic effect. We were able 



to measure it for both filters and lenses of the ISO-HDF 
observation. 

The flat fields are computed in CIA by making the 
assumption that the illumination of the detector is ho- 
mogenous. Because of the field distortion, this is not truly 
the case since the pixel size is not uniform on the sky (the 
pixels at the edges of the array cover a wider surface). 
A new flat-field correction has to be applied in order to 
account for it. This flat-field is of the form : 



S 



16,16 



where Sij is the surface on the sky of pixel i,j. The 
pixel (16,16), which is the center pixel of ISOCAM LW 
array and therefore the least distorted, has been taken as 
reference. The images should therefore be divided by F 
before reprojecting on the sky. 



4--4- Astrometric corrections 

As the observations of the ISO-HDF have been performed 
in three separate rasters, one has to coadd all of the images 
in order to reach the maximum sensitivity. But this can 
only be done if the pointing of each raster is well known. 
This is not straightforward because of : 

1. The absolute pointing accuracy (i.e. the 1 a error that 
can be made when acquiring a new target) of ISO is of 
3 arcsec. This is much greater than the relative point- 
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ing accuracy, (i.e. the accuracy of pointings between 
two raster steps), which is less than 0.1 arcsec. 
2. The lens wheel jitter. In order to avoid any mechan- 
ical blocking, the gear wheel has been designed with 
a small play. The position on which the lens stops is 
not fixed. It has been shown by the Cam Instrument 
Dedicated Team (CIDT) and by us that there are only 
two positions that each lens can take for a commanded 
position, probably at either side of the play. This can 
be very easily detected by a close inspection of the flat 
field derived from the data. If the leftmost column of 
the detector receives very little light, we know we are 
in the 'left' position. This jitter results in an offset of 
about 1.2 pixels from the optical axis, thus « 7" with 
the 6" PFOV. It modifies also the distortion pattern of 
ISOCAM that has been measured for both positions. 
The only way to avoid such problem when coadding 
various rasters is to use the screw 129 and have them 
observed in a sequence in order to prevent any move- 
ment of the lens wheel between the two rasters. This 
was not the case for the ISO-HDF observations. 

If all the observations suffer from the absolute pointing 
error, only one raster of the LW3 observation suffers from 
the second problem. The only way to coadd efficiently 
the rasters is to identify common sources in each of them 
and to compute an offset. These sources have to be bright 
enough to have a well definite position. In order to correct 
the absolute pointing error, good astrometric references 
are needed. The astrometry of the HST observation of 
the HDF has been established very carefully (Williams et 
al. 1996), but is not useful in our case : since the FWHM 
of ISOCAM PSF is large, there are in general several op- 
tical counterparts which could be associated to a single 
ISO-HDF source. 

In order to solve this problem, we used the well known 



Table 2. ISO-HDF astrometric correction. The radio 



sources (Fomalont et al. 1997) used are listed for each 
field. 



MIR/radio correlation (Helou et al. |1985| ). Since the HDF 
and the flanking field s have been observed at 8.4 GHz 
(Fomalont et al. 1997 ), we used this source list to match 
our brightest objects. We were successful in identifying at 
least 4 radio sources in each LW3 raster and at least one 
in each LW 2 raster. These sources are given in Table |] 
together with the derived offsets. 

4-5. Detection, source photometry and transients 



As explained in Starck et al. ( |1998| ), the source detection 
is made on the final reprojected map, against the noise 
map, using a wavelet technique. We use a b-spline wavelet 
transform, and detect at each scale all pixels above a given 
threshold (t w ) of the transformed noise map. In the ISO- 
HDF case, the final size of the pixels is 3" in LW3 (1.5" in 
LW2) for a PSF that is 9" FWHM (resp. 4.5). We therefore 
ask the algorithm to perform the detection down to the 
3 rd scale, that is over a box of 16 pixels, matching the 
PSF size. Since the maps are not square, some sources are 



Raster 


8.4 GHz 


a offset 


S offset 


LW2 Raster 2 


3646+1404 
3649+1313 


4.5" ± 2.0 


+3.0" ± 2.0 


LW2 Raster 3 


3649+1313 


0.0" ± 2.0 


0.0" ± 2.0 


LW2 Raster 4 


3644+1133 


+4.5 "± 2.0 


+2.5" ± 2.0 


LW3 Raster 2 


3634+1212 
3634+1240 
3644+1249 
3646+1404 
3646+1448 
3649+1313 


0.0" ± 0.5 


+4.0 " ± 0.5 


LW3 Raster 3 


3644+1249 
3646+1404 
3639+1313 
3653+1139 


-0.7" ±0.5 


+4.0 " ± 0.5 


LW3 Raster 4 


3634+1212 
3634+1240 
3644+1249 
3646+1404 
3646+1448 
3649+1313 
3653+1139 


0.0" ± 0.5 


-4.2 " ± 0.5 



detected at the edges of the field due to non continuous 
limit conditions. Such sources are eliminated. 

It is very important to note that the t w that is em- 
ployed here describes the noise level of the detection map, 
but that this noise is not Gaussian. Since many treatments 
applied on the cube where done using median transform, 
and since the 'dippers' and 'faders' are highly non Gaus- 
sian, this means that one cannot relate a3r m detection to 
a 99.69 % confidence level, as it is the case for Gaussian 
statistics. Simulations only can give confidence levels for 
the value of the detection threshold. 

The photometry of the objects can be measured in two 
ways. First, the detection program constructs an image of 
each object. The total flux in this image is computed. 
Second, a direct aperture photometry can be performed 
on the ISOCAM map. We have used the second method 
because simulations have shown us that for faint objects, 
the reconstruction using wavelets is not reliable. The mea- 
sured Analog to Digital Units (ADU) are corrected to 
account for transients and for PSF. This correction has 
been obtained by simulations (see below). The corrected 
ADU are then translated to milli Jansky by using ISO- 
CAM 'cookbook' conversion tables. 

The response of the ISOCAM LW detector is affected 
by a transient when its illumination changes. The response 
reaches instantaneously a level corresponding to 60 % of 
the final, stabilized signal, the rise for the remaining 40 % 
being a slow exponential whose time constant depends 
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on the past history of the pixel. First order corrections 
have been produced, especially by Abergel et al. ( 1990| ). 
Transients are the source of two main problems : 

1. an error in the photometry if the measured flux is not 
stabilized. 

2. 'ghost' sources appear on maps after that the detector 
has observed a bright source (> 500/iJy). 

Since the ISO-HDF has very small raster steps and no 
bright source is present in the field, the second point is not 
a problem in our case. Correcting transients has an effect 
on the time history of each pixel. Some glitches, with their 
slow decrease, mimic a transient behaviour, and could be 
treated inadequately by the correction algorithm. Since we 
detect these effects by criteria on their shape in the wavelet 
space, we might miss some of them after the correction of 
the transients. On the other hand, transient effects can be 
modeled, and taken into account when converting ADU to 
physical units. For these reasons, we do not perform any 
transient correction during the cleaning of the data cube. 

Simulations were widely used to calibrate the transla- 
tion of measured ADUs in 'real' ADUs, then in fj,Jy and 
arc detailed below. 



5. Simulations 

The PRETI method has been developed to detect faint 
sources in all ISOCAM data (if they are taken in fields 
where the background is more or less constant), but has 
to be tuned in order to work very efficiently, as shown 
in the previous section. This fine tuning has been done 
through various tests and simulations, that have proven 
to be very useful to ascertain detections, study the com- 
pleteness, the number of false detections and to calibrate 
the source photometry. 

In order to test the PRETI method, we have tried to 
build simulated data. These should be as near as possible 
of real ISOCAM data, that is, contain glitches, transients, 
and sources. A special attention has to be payed to the 
'fader' and 'dipper' long term effects of glitches since they 
are the most limiting factor for the faint source detection. 
Unfortunately, no available model is able to reproduce ac- 
curately these behaviours of the LW detector. The only 
solution is to use real ISOCAM data to build simulations. 



5.1. cosmic rays in simulated data 

Cosmic ray effects are dependent of three parameters 

1. the integration time because it is directly related to the 
mean number of cosmic rays that will hit the detector 
per readout. 

2. the background emission, because in the long term 
effect a cosmic ray may affect the gain of the pixels 
(though this is not yet clearly established). 



3. the period of the observation since the cosmic ray den- 
sity varies along the orbit of the satellite as it crosses 
the radiation belts. We have noticed variations of 40 % 
in the number of 'faders' and 'dippers' between vari- 
ous observations, early part of the orbit being most 
affected. 

For these reasons, data for the simulation are taken 
with the same integration time, same gain, same filter and 
zodiacal background as the observation that is to be sim- 
ulated. Moreover, in order to avoid any modification of 
the time history of pixels, these data have to be taken on 
periods at least as long as the observation. For example, 
data taken on a shorter period and duplicated would cause 
discontinuities in the time history and produce spurious 



5.2. Suppression of real sources in the data set used for 
simulations. 

Since we look for faint sources, it is very likely that the ob- 
servations selected to be used for simulations will contain 
such sources. 

The ideal case is when these are 'staring' observations, 
i.e. when the satellite points during a very long time at 
the same position. Since the flat-field is built directly from 
the observation, any source present in the field is cleaned 
during the operation. 

One less ideal case is when the raster parameters of 
the observation are different from the simulated data, the 
sources in the simulation data set are blurred when the 
raster map is reconstructed with different steps. If there 
are bright sources in the simulated data, which would de- 
spite the blurring process still have a high S /N level, these 
can be masked out. 



5.3. addition of sources 

Point sources of various fluxes are added using the PSF 
model of M. Perault and K. Okumura, that match well the 



measurements (Okumura et al. 1997). Each source flux 
is converted in ADU by using the ISOCAM cookbook. 
Object are randomly distributed on the raster map, at 
a precision of a tenth of a pixel, then each position on 
the detector array is recomputed, without taking the field 
distortion into account. 

Each raster image is then translated into a data cube 
by using the transient model from Abergel et al. (1996), 
where the detector response I? is a function of the incom- 
ing flux / with : 



D(t) = rl(t) + (l-r) / /(t')exp 



t-t 
a/ '1(f) 



dt' 



It is clear from this formula that the response depends on 
the intensity of the source and on the background, because 
the exponential time constant is a function of I. For this 
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Input flux ( M Jy) 



(a) LW2, 5t w detection level. 




Input flux (fUy) 

(b) LW3. Solid : 5t w detection level. Dashed : 7t w detection 
level. 



Fig. 6. Percentage of detected sources in simulations of the ISO-HDF, as a function of the input flux of the sources, 
for both filters 



reason, we add to the simulated source images the me- 
dian of the deglitched cube that is to be simulated before 
running the transient simulator. We subtract this image 
later, in order to have a simulated cube on a zero back- 
ground, and add it to the raw data selected to perform 
the simulation. 

5.4- ISO-HDF simulations 

We have searched the ISOCAM calibration database for 
observations matching the ISO-HDF configurations. For 
the LW3 observation, we found a set of flat-field obser- 
vations, called 'FLAT16' that were done with the same 
filter, the same lens and same integration time. Seven 
bright sources present in the field were masked out. Fainter 
sources have been blurred since this observation is a raster 
with half an array of raster step (16 pixels). We checked 
the PRETI method by running it completely on this set. 
We set our deglitching parameters in order to avoid any 
detection in that field (since their should be no source), 
except the ones that could be related to a source present 
in the original dataset. We also checked that the 'fader' 
and 'dipper' rate is about the same in both observations. 
This set of data is a perfect tool to build complete sim- 
ulations. Since it is an empty field, we are able to study 
very carefully the level at which PRETI begins to give us 
false detections. 

For the LW2 ISO-HDF, we could not find any sat- 
isfactory set of data : only a few calibration have been 
performed with an integration time of 10 s, and they are 
always short ones. We eliminated also data from our ID- 
SPCO survey program that were matching all the con- 
figuration requirements, because we found a rate of long 



term effect cosmic rays about 20 % higher : this causes 
severe differences in the sensitivity that can be reached 
by an observation. Since our first reductions show us that 
the source density is not very high in the ISO-HDF LW2 
observation, we decided to use these data themselves for 
the simulations. The drawback is that this prevents us to 
fully study the confidence level and the false detections. 

We performed two kinds of simulations. First, 10 
sources of a given flux were randomly added, and the op- 
eration was repeated from 1 mJy to 100 /xJy by steps of 
100 /iJy, and from 100 to 10 fiJy by steps of 10 /xJy. This 
allowed us to investigate our detection limits and to cali- 
brate our source photometry. Finally, we performed sim- 
ulations containing a wide range of fluxes, more or less 
matching number count models, in order to test the com- 
pleteness of the method. 

5.5. Simulation results on the ISO-HDF 

Figure ^ presents the results of our counts on simulated 
data for LW2 (a) and LW3 (b). For the LW2 filter, all 
added sources in the field are detected down to 60 //Jy 
for a 5 t k detection, showing that we are complete at 
this level. However, this result does not take into account 
confusion. A faint source is not detected if is is too close 
to a bright one ( ~ 9" in LW3, depending on the relative 
brightness of the two sources) . This effect has to be taken 
into account when deriving number counts. For the LW3 
observation, the 99.9% completion is obtained above 200 
/zJy, for a 5t w detection and stays above 85% at 100 fiJy. 
7t w detection, that are more reliable in terms of avoiding 
false detections, are only 99.9% complete above 1 mJy, but 
stay of the order of 90 % down to 220 /iJy. We could not 
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(a) Measured flux as a function of input flux for 2 photome- 
try methods. Solid : 2 pixel aperture photometry (3"). Dot- 
dashed : wavelet photometry. Dashed : the identity curve. 



(b) Ratio of measured flux over real flux as a func- 
tion of the real flux. Solid : 2 pixel aperture photome- 
try (3"). Dot-dashed : 3 pixel aperture photometry (4.5"). 
Dashed : wavelet photometry. 
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(c) Measured flux as a function of input flux for the wavelet 
method. Solid : measure. Dash-dotted : la error on mea- 
sure. Dotted : Identity. 



(d) Measured flux as a function of input flux for 2 pixel (3") 
aperture photometry. Solid : measure. Dash-dotted : la 
error on measure. Dotted : Identity. 



Fig. 7. Photometry results of LW2 simulations. Comparison of photometric methods and errors. 



measure the number of false detections in the LW2 band 
since we used the data set itself to do the simulations. 
Since PRETI is run with the same arguments on both sets, 
we assume that the difference should be small between the 
two filters and tune our criteria according to LW3. 

Figures ^ and || show the calibrations curves deduced 
from simulations, i.e. the relation between the measured 
flux to the real flux that was simulated, for two mea- 
surement methods : (1) an aperture photometry taking 
a radius of 3" for LW2 and 6" for LW3 (two pixels in 
each case), (2) the integrated flux of the objects recon- 



struct ed by our wavelet based detection program (Starck 
et al. 1998 ). The two methods display very different be- 
haviors for the faint end of the plots. 

Photometry on wavelet reconstructed objects (panels 
(c)) shows a non-linearity effect for input fluxes below 100 
/iJy in LW2 (resp. 200 /zJy in LW3). Above these lev- 
els, the ratio of the measured flux over the input flux is 
roughly constant (~ 0.45 for LW2). We interpret this fac- 
tor in terms of the effect of transients (factor of 0.6) , 
and by the fact that only the central core of the PSF is 
reconstructed by the algorithm (factor of 0.8), thus giv- 
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(a) Measured flux as a function of input flux for 2 photom- 
etry methods. Solid : 2 pixels aperture photometry (3"). 
Dot-dashed : wavelet photometry. Dashed : the identity 



(b) Ratio of measured flux over real flux as a function of 
the real flux. Solid : 2 pixel aperture photometry (3"). Dot- 
dashed : wavelet photometry. 
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(c) Measured flux as a function of input flux for the wavelet 
method. Solid : measure. Dash-dotted : la error on mea- 
sure. Dotted : Identity. 



(d) Measured flux as a function of input flux for 2 pix- 
els (3") aperture photometry. Solid : measure. Dash- 
dotted : la error on measure. Dotted : Identity. 



Fig. 8. Photometry results of LW3 simulations. Comparison of photometric methods and errors. 



ing a theoretical correction of 0.48. Below these levels, 
the non-linearity is explained by the fact that the wavelet 
analysis program reconstructs objects with extended low 
levels wings. This effect might be due to the fact that we 
use a b-spline wavelet transform of our map that does not 
match closely enough the ISOCAM PSF. We note that the 
non-linearity is more marked for LW3 than for LW2 where 
the PSF is less extended, and is therefore closer to the 
wavelet shape. For the last filter, the measured level stays 
constant at 30 /iJy while it rises in the first one for de- 
creasing fluxes. The reason for this difference is not clear, 



but might be explained by the confusion of the the ob- 
servations in the LW3 map, or by residual glitches^]. This 
non-linearity prevents to use the "wavelet" photometry 
for low levels objects, thus this method is not appropriate 
for the HDF data set. 

Aperture photometry (panels (d)) does not display 
such non-linearity. In LW2 , the ratio of the measured 
flux over the input flux decreases as the input flux de- 
creases. The explanation for this behaviour is the follow- 



1 the only way to test this hypothesis is to achieve clean 
simulations of empty fields ! 




Fig. 9. Map of the ISO-HDF LW2 filter (6.75 /xm). Resolution is of 1.5" /pixel. Contours are 5 t w detections levels at 
35, 50 and 100/jJy 



ing : when the flux decreases, the outer pixels of the cen- 
tral core of the PSF are more and more dominated by 
the noise, and the summed flux contains therefore nega- 
tive terms. Indeed, due to the baseline substraction, a few 
readouts are slightly negative before and after the source, 
(the background level is overestimated). Therefore, the av- 
erage level on the final map around the source is slightly 
negative. This effect is illustrated in figure |7] (b) were a 
3 pixel radius aperture photometry has also been drawn 
(dash-dotted curve). For this method, the ratio at high 
flux is larger than for a 2 pixel radius aperture (more of 
the central core is taken into account), but decreases more 
abruptly toward faint fluxes, because more noise and neg- 
ative values are integrated in the aperture. We note that 
at high flux, the ratio between the measured flux in the 3 
pixels and 2 pixels apertures is ~ 0.3/0.25 = 1.2, in good 
agreement with the PSF model that predicts 1.209. Again 
for this photometric method, the behaviour of the LW3 fil- 
ter is slightly different, because the slope of the calibration 
curve is less steep toward faint fluxes. The change of slope 



occurs at the same level (200 /xJy) where the "wavelet" 
method begin to be non-linear. This indicates that the 
same process is at work for both methods, but that aper- 
ture photometry is less affected. For these reasons, we have 
used a two pixel aperture photometry for the HDF. 

6. Results on the ISO-HDF 

Figures [)] and [n] present the reduced map of the ISO- 
HDF LW2 and LW3. Many sources are clearly visible on 
the latter, while the former appears almost empty. Var- 
ious sensitivity levels have been overlayed on the maps, 
showing that the noise is very inhomogeneous. These sen- 
sitivity levels were computed by comparing the noise level 
at each scale with the wavelet coefficients of a point source 
image, the source being centered on a pixel. Thus, these 
levels give only an indication. If a source fell between two 
pixels, its wavelet coefficients at lower scales (high spatial 
frequencies) is lower, while those at higher scales (low spa- 
tial frequencies) are higher. Since the noise decreases when 
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Fig. 10. Map of the ISO-HDF LW3 filter (15 fim). Resolution is of 3" /pixel. Contours are 5 t u 
35, 50, 100 and 200/xJy 



detections levels at 



considering higher scales (Starck et al. 1998), sources can 
in fact be detected at fluxes lower that the levels shown 
on figures [)] and [l(]. 

6.1. Source catalog 

Table | gives the catalog of detected sources at 6.7 /im 
(LW2 filter) and 15 /xm (LW3) above a secure 7t w of pho- 
ton and readout noise. The first column is our ISO-HDF 
identificator, the second and third column are the right 
ascension and declination of the object. The pointing ac- 
curacy is of 1 pixel on the final map, i.e. 3" for the LW3 
image and 1.5" for the LW2 image . Columns 4 and 5 give 
the objects fluxes in the LW2 band in fj,3y and ADU/G/s, 
that is in analog to digital unit per second divided by the 
instrument gain, and columns 7 and 8 give the same in 
the LW3 band. Columns 6 and 9 give the detection level 
in the LW2 (resp. LW3) band. Column 10 identifies the 
field where the object is observed (HDF or flanking field). 
Column 11 gives the spectrometric redshift when avail- 



able. Finally, column 12 give the morphological type of 
the object. 

The redshift and morphological type of each ISOCAM 
source depend on the identification of its optical counter- 
part. Hopefully, our astrometric correction is good enough 
to allow us a good superposition with optical images, at a 
scale of 1 pixel. This is clearly visible on figures 11 and [l2] 
where we have overlayed our ISOCAM maps on the KPNO 
IRIM deep K exposure ^ (Dickinson at al. in prep), as well 
as on figure |l3| where we have overlayed the LW3 ISOCAM 
map and the LW2 detections on the color composite image 
of the HDF observed by HST from williams et al. (|l996fl. 



In most cases, the ISOCAM source matches an IR or op- 
tical source, and the identification is easy. However, due 



2 Based on observations made at the Kitt Peak National Ob- 
servatory, National Optical Astronomy Observatories, which is 
operated by the Association of Universities for Research in As- 
tronomy, Inc. (AURA) under cooperative agreement with the 
National Science Foundation. 
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Table 3. PRETI main source list for the ISOHDF observations at 6.5/zm (LW2) and 15/xm (LW3). All sources are 
detected at 7t w above readout and photon noise. (1) : Source identification 'PM3' is for PRETI Main list LW3, 'PM2' 
for PRETI Main list LW2 . (2), (3) : source position for the J2000.0 equinox. (4) : source flux in micro Jansky in the 
LW2 filter with error bar or upper limit (5r w ). A '-' indicates that the source is not in the LW2 field. (5) : source flux 
in ADU/G/s in a 3" aperture. (6) : detection level (nr w ) in the LW2 image. (7) : source flux in microJansky in the LW3 
filter with error bar or upper limit (5t w ). (8) : source flux in ADU/G/s in a 6" aperture. (9) : detection level (nr w ) in 
the LW3 image. (10) : Identificator of the HST field in which the source falls (hd stands for HDF, others for flanking 
fields). (II) : spectroscopic redshift with source indicated by a note, (c) : data from Cohen et al. ( |1996| ), (p) : data 
from Phillips et al.( 1997| ), (i) : compiled in Cowie et al. (1997), ( h) dat a from Hogg et al. (1977b). (12 ) : mo rphological 
type of the object, '(i)' indicates a type compiled in Cowie et al. ( 1997| ),'(p)' type from Phillips et al. ( |l997|) , otherwise 
determinated by us. E : Elliptical or SO galaxy. S : spiral. M : merger. P : peculiar/irregular, C : compact object from 
Phillips et al. ( 1997|) sample, A : active nuclei galaxy, G : source falls on a group. U : unknown type. N : no optical 
counterpart. A '?' indicate a dubious type identification. 
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Fig. 12. Comparison near-IR/LW3 filter. The contours of 
the central part of the LW3 map have been overlayed on 
the IRIM K map. Contour levels begin at 1.5 /iJy with 
steps of 2 /xJy. 



to the large PSF (15" at 15/zm) of ISOCAM, sometimes 
more than one object could be associated with our detec- 
tion. Figure |TJ displays an example of such a case : the 
main peak of source HDF_PM3_6 correspond to two ob- 
jects. In this particular case, we have taken the brightest 
object in K as counterpart. When K data are not avail- 
able, we have used the 'C (for group) symbol in our table, 
because we could not determine the right counterpart. 

Looking closely at figure |l4|, one can see that the object 
HDF_PM3_6, detected as a whole, displays substructures. 
Our detection program does not separate the various com- 
ponents of the source because it is not robust against 
blending : to avoid detection of cosmic ray residuals, we 
choose to perform the detection starting from the second 
scale of the wavelet transform of our map, i.e. on scales 
of the order of 3 pixels that prevent us to separate very 
closeby sources, especially if they are much weaker than 
the main detection. However, in some cases, this separa- 
tion can be done by hand, when blending is weak enough 
to allow for the separation of the source (Rayleigh crite- 



ria). This is the case for source HDF_PM3_6, for which 4 
extensions are visible on figure [l4|, South and East of the 
main peak. These extensions are in our supplementary list 
(see table g) with the names HDF_PS3_6a, HDF_PS3_6b, 
etc.... Figure |lj displays another case where the sepa- 
ration is possible, for source HDF_PM3_17 . Figure [l(] 
displays a closeup on source HDF_PM3_27 : the source 
is extended, with multiple counterparts, but we failed to 
separate the various contributions. We note however that 
this source is the radio source 3639+1313 in the catalog of 
Fomalont at al. (1997), whose contours in the radio map 



(Richards et al. 1997) are also clearly extended. 

We produced a supplementary list of sources with these 
extensions of reliable sources, together with results of 
PRETI at lower detection levels (6 and 5 t w ), that is listed 
in table |j. However, this catalog may contain spurious de- 
tections, as it is certainly the case for sources HDF_PS3_13 
and HDF_PS3_29, with fluxes below 10 /iJy. In the future, 
we expect to obtain data of a long staring observation on 
an empty field. Using this observation as a basis for simu- 
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Fig. 13. Overlay of the ISO-HDF observation on the HDF color composite image produced by Williams et al. (1996 
)from the F450W, F606W and F814W images obtained with HST. Orange contours : contours of the 15 [im (LW3) 
ISOCAM observation. Contours levels begin at 1.0 /iJy with steps of 2 /zJy. Green circles : ISOCAM 6.75 /im (LW2) 
detections. 



lations, we hope to be able to obtain higher reliability on 
the sources at these low level of flux. 

6.2. Catalog properties. 

Our main source list consist in 49 objects detected above 
a 7t w detection threshold in both filters. 42 sources are 
detected at 15 (urn only, 4 at 6.5 and 15 /im, and 3 at 



6.5 /im only. Since the LW3 field is larger than the LW2 
field, this last number means that 21 of the LW3 sources 
are not detected in LW2. The supplementary list adds 
51 sources to this catalog for a total of 100 objects. 47 
are observed in the LW3 filter only, and 4 in the LW2 
filter only. None are detected in both filters, 15 of the 
LW3 sources that could be visible in the LW2 field remain 
undetected. For these missed detection, we have computed 
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HDF_PM_17. The contours of the details of the LW3 map 
around source 17 have been overlayed on the IRIM K map. 
Contour levels begin at 1.5 fiJy with steps of 2 fiJy. 



5t w upper limit as follows : a test image containing only 
the source is computed by projecting a point source as 
observed at each raster pointing in a copy of the final 
raster map. This ensures that the modification to the PSF 
produced by microscan are taken into account. The b- 
splinc wavelet transform of this image is then calculated 
and compared to 5 times the noise map at each scale. The 
minimum flux ensuring the detection is then converted in 
microJansky, using our calibration from simulations. 

All sources from the main list within the HDF have op- 
tical counterparts, thanks to the depth of the HST map, 
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Fig. 16. Comparison nea r-IR/LW3 filt er for source 27. 

The contours of the details of the LW3 map around source 
HDF_PM_27 have been overlayed on the IRIM K map. 
Contour levels begin at 1.5 [iJy with steps of 2 /iJy. 



but 6 sources within the flanking fields do not have vis- 
ible counterparts in the much less deep F814W images. 
However, these 6 sources are very reliable and were also 
found by other teams (see below) . The number of sources 
without counterpart rises to 11 in the supplementary list, 
of which 1 lays in the HDF. Source HDF_PS3_6e does not 
have any counterpart in the 'iw' HST image, but has a 
counterpart in the IRIM K image. This example shows 
that the reliability of the source should not be established 
on the existence of counterparts in optical images only, 
especially in the flanking fields. It is remarkable that all 
but one (3652+1444) of the radio sources observed at 8.4 
GHz by Fomalont et al. (1997) in the field of view, are 
detected, mostly in the LW3 filter (Table g), and one in 
the LW2 filter only (Table g). 

Other teams have used different algorithms to produce 
source lists. We have compared our results to those of 
the Imperial College group (hereafter IC) (Goldschmidt 

|l99|). 



The 



et al. 1997) and to those of IAS (Desert et al. 
common source lists are given in Table g for the LW3 filter, 
and in Table g for the LW2 one , as well as the names of 
the radio sources observed at 8.4 GHz by Fomalont et al. 
( |l99"7j ) matching ISOCAM sources. 

PRETI a nd the "triple beam switch method" of IAS 
(Desert et al. 1998) give very similar results (see Table g 
and g), especially in the LW3 band where out of the 41 
sources detected by them, we find 37. We do not detect 
HDF_ALL_LW3_8, 14,29,40 that are at 3.4cr, 3.3cr, 3.9cr 
and 3.1cr respectively, that is among the faintest of their 
list. Both number counts agrees above 200 /xJy where 
PRETI reaches 99.9% completeness : the thwo methods 
detect 24 sources. Figure |l7| displays a comparison be- 
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tween the fluxes derived by the two methods for the com- 
mon source list. The fluxes derived by the "triple beam 
switch method" are lower than our fluxes by a systematic 
factor of 0.82. We note that both methods relie on the 
ISOCAM cookbook values for the correction of ADUs into 
milli-Janskies. For one part, we can explain the difference 
by an implicit color correction that we have introduced 
in our simulations : our PSFs are simulated assuming a 
rising spectrum with spectral index o f 3, w hile standard 
ISOCAM PSFs, used by Desert et al. 



1998 



are measured 

on stars, thus assuming a spectral index of (-2). When 
comparing our PSFs to standard ones, a mean difference 
of 10 % is found. Moreover, source fluxes are measured 
by us with aperture photometry while the "triple beam 
switch method" uses fits of gaussian of fixed width. Since 
at 6", the ISOCAM PSF is only roughly approximated 
by a gaussian, this can account for the remaining discrep- 
ancy. At these faint level of fluxes, such differences are 
not very important, considering that the source lists of 
the two catalogues matche very well. However, PRETI al- 
lows to reach deeper sensitivity, as we consider sources at 
40 ^tJy as reliable. For the LW2 filter, the "triple beam 
switch method" detects 6 sources, 3 of which are also de- 
tected with PRETI. Of the 3 undetected objects, 2 lay on 
the edges of the map where we do not allow for detection. 
The last one is measured by Desert et al. at 58 /iJy where 
we are not complete. 

The IC group have made the original data reduction of 
the ISO-HDF, when some of the latest refinements of ISO- 
CAM calibration such as field distortion were not avail- 
able. We find 21 of the 22 sources detected by them in 



their LW3 image (Goldschmidt et al. 1997). We could not 



find the source J123659. 4+621337 in our map. Source pho- 
tometry differs from our catalog : on the average, their 
flux estimates are higher than ours by a factor of 1.5, as 
shown in figure O. They performed simulations, and con- 
cluded that they should be 70% complete at a level of 
225 /iJy. According to the 1.5 factor of difference in pho- 
tometry, they should have detected at least 26 of our 37 
sources above 150 [iJy, where they only find 19. This indi- 
cates that they are only 50% complete at this level, while 
PRETI reaches 95% at a level of 150 /zJy, according to 
our simulations, if source confusion is ignored. Comparing 
the LW2 detections, we only find 5 of the 27 sources in 
their list (3 from their main list of 6 objects and 2 from 
their supplementary list of 21 objects), as shown in ta- 
ble ||. Moreover, we measure much fainter fluxes for these 
objects. Since this observation is made with 10s integra- 
tion time, so that pixels affected by cosmic ray glitches are 
more numerous, two explanations can account for this : ei- 
ther we set too heavy a deglitching criterion and we have 
cleaned real sources out, or they did not clean their data 
enough against cosmic rays. Our analysis lead us to favor 
the second explanation, because we do not find that we 
are missing sources in the simulations. Here also the re- 
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Fig. 17. Comparison of the photometry of common 
sources in the ISO-HDF at 15 /im for PRETI and IAS. 
Solid line is 1 to 1 relation, dashed lines are a 20 % error 
with respect to the 1 to 1 relation. 



suits of Desert at al. (199S) are in good agreement with 
us. 

6.3. Number counts 

The small number of LW2 sources does not allow to com- 
pute reliable number counts, as the poissonian fluctuations 
introduce an uncertainty of plus or minus 40 % 1 a when 
using the reliable 7t w source list. Nevertheless, we have, 
at the level of 100 /iJy, 5 reliable detections in an area of 
10.45' 2 , that is : N(> lOO^Jy) = 1.7 x 10 3 ±8 x 10 2 sources 
per square degree, a value compatible with the predic- 
tion of a pure passive lumi nous e volution model of number 
counts (Franceschini et al. 1997 ), as shown on Fig. [19[ 

The more numerous detections in the LW3 band allow 
us to compute some points of the Log(N)-Log(S) curve, 
and to try do determine its slope. We determined two 
curves, using the 7t w and 5t w detection thresholds, for 
which a few more spurious detections are expected. 

In order to compute reliable number counts, we have 
taken into account the following parameters : 

1. the area available for detection at a given flux with a 
given detection threshold (Nt w ) has been computed, 
using the same technique as for determining the upper 
limits of undetected sources. However, an additional 
hypothesis was made : each source is considered cen- 
tered on a pixel of the map. 
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Fig. 18. Comparison of the photometry of common 
sources in the ISO-HDF at 15 for PRETI and IC. 
Solid line is 1 to 1 relation, dashed lines are a 20 % error 
with respect to the 1 to 1 relation 

2. the blending of faint sources by bright ones is treated 
in a very simple way : the area occupied by the sources 
of all the brighter bins of flux is subtracted from the 
area available for the detection of sources at a given 
bin. For each source, we have taken out a circle of 
radius 9" , that is approximately the distance between 
sources separated by PRETI. 

3. we take the completeness of the detection at a given 
level into account by using the results of our simula- 
tions (see Fig. ^j). 

4. when dealing with the 5t w detections, we have not 
taken source HDF_PS3_40 into account, because of its 
high probability of being a star. 

This method leads us to derive the number counts that 
are plotted on Figs. ^ and gl], with error bars correspond- 
ing to lcr poissonian fluctuations (as the square root of the 
number of events). Moreover, we have also plotted on these 
figure two extreme limits : 

1. The 'upper' limit, where all sources are given their 
measured flux plus their positive error bar, to which 
we add the poissonian fluctuations. 

2. The 'lower' lower, where all sources are given their 
measured flux minus their negative error bar, to which 
we subtract the poissonian fluctuations. 

Finally, we have also plotted on these figures the pre- 
diction of a model with no evolution (Franceschini et 
al. |1997|) . 



Fig. 19. Number counts in the ISO-HDF at 6.75 /Ltm, 

based on the 7t w sample. 

Solid : no evolution model from Franceschini et al. (1997) 
Hatched area : lcr error area for the integrated counts 
point (cross) from the 7t w sample. 



In both 7t w and 5t w sample, the straightforward result 
is that the counts derived from ISO-HDF in the LW3 band 
show a strong excess : nearly a factor of ten at 200 /iJy, 
with respect to the prediction of the no evolution model. 
This e xcess is higher than the one obtained by Oliver et al. 
( |1997| ) from the first reduction of the ISO-HDF data. This 
might be due to the underestimation of the completeness 
of the sample. 

In order to test our simple method to account for 
blending, we have plotted in figure ^l] the number counts 
obtained with the 5t w sample with the method described 
when blending is not taken into account, but adding to our 
sample the unblended sources in our supplementary list. 
On this last plot, the slope of counts is constant down to 
60 /xJy, and slightly steeper than predicted in the frame- 
work of the no evolution model. 



7. Discussion 

The number counts derived from the ISO-HDF observa- 
tions implying strong evolution at 15 /im and no evolution 
at 6.75/im, may appear surprising at first sight. This ap- 
parent discrepancy is in fact easy to interpret, considering 
that we are not observing the same physical processes at 
these two wavelengths. 

As shown in Fig. ^, the LW2 filter tends to select the 
old stellar population of elliptical galaxies; taking advan- 
tage of the K correction, the detection of such systems 
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Fig. 20. Number counts in the ISO-HDF at 15 /im, based 
on the 7t w sample. 

Solid : 7t w with poissonian error bars, 
dashed : 'upper' and 'lower' limits of counts, 
dash-dotted : no evolution model 



are favored at z > 0.5. For gaseous systems, dominated by 
dust and UIB emission, as the redshift rises, the LW2 filter 
crosses a gap between the UIB features above 6.2 /jm and 
stellar emission below 4/mi in the rest frame, leading to a 
small number of detections, even though the sensitivity is 
higher than with the LW3 filter. 

The LW3 filter picks up hot dust in nearby systems as 
well as redshifted UIB emission at z as high a 1.5. Studies 
of nearby galaxies show that this emission is enhanced in 
the regions surrounding a starburst, so that the detection 
of starbursts, often in mergers, is favoured. This is borne 
out in the results. Indeed, a close examination of the mor- 
phology of ISO sources shows that early type galaxies are 
numerous among the few LW2 detections (6 out 11 de- 
tections), while spirals, irregulars, or mergers dominate in 
the LW3 sample. 

The z distribution of our sample of LW3 sources (with 
signal in wavespace 7 times above the threshold) is dis- 
played in Fig. |2^, and compared to that of sources with 
B < 25 and K < 22., taken from the list by Cowie et 
al. ( |1997| ). The redshifts in our sample range from z = 
0.078 to z=1.242, with a median value of z=0.585 ; most 
of ISOCAM LW3 sources have z > 0.4 . The distribution 
in redshift of our two samples is very similar, implying 
that the LW3 filter does not select a particular range in z. 

A forthcoming paper (Elbaz et al., in prep.) will 
present a detailed study of each ISO-HDF sources. This 
study will allow to determine if the evolution observed at 
15 /im is a general trend in galaxies or is due to certain 
types in the galaxy population, and put limits on the evo- 
lution of the early type systems at z <1.5. 




Fig. 21. Number counts in the ISO-HDF at 15 /im, based 
on the 5t w sample plus unblended sources. 
Solid : 5t w with poissonian error bars, 
dashed : 'upper' and 'lower' limits of counts, 
dash-dotted : no evolution model 



8. Conclusion 

The PRETI method for faint source detection with ISO- 
CAM has been successfully applied to the observation of 
the HDF and flanking fields. Faint fluxes are reached with 
a good reliability. Moreover, our simulations allow us to 
test the completeness and photometric accuracy of our 
results. On the ISO-HDF, the method gives a 99.9% com- 
pleteness at 200 /xJy in LW3 when detecting above five 
times the noise level in wavelet space. The same com- 
pleteness is obtained at 65 /xJy in LW2. Using PRETI, 
we detected 49 sources in our main source list, plus 51 
in a less secure supplementary list. Most of the sources 
are detected at 15 /im, 11 are detected at 6.7 /im only. 
Our results for the brightest sources are in good agree- 



ment with those of Desert et al. (1998) in both filters, as 
well as with those of Rowan- Robinson et al. (1997) at 15 
/mi, but not at 6.75/^m. The number counts derived for 
the LW3 filter show a clear excess of sources with respect 
to the predictions of a no evolution model (Franceschini 



et al. 1997). A forthcoming paper (Elbaz et al. in prep) 



will refine these results with a source by source study and 
discuss extinction and star formation rates. 
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Fig. 22. Histogram of redshifts distribution in two sam- 
ples 

Solid : ISOCAM LW3 sources . 



dashed : Cowie et al. (J1997J) sample on HDF with B < 25 
and K < 22. 
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Table 4. PRETI supplementary source list for the ISOHDF observations at 6.5/Ltm (LW2) and 15/im (LW3). All LW3 
sources are at least detected at 3^ above readout and photon noise, while all LW2 sources are detected atleast above 
5t w . (1) : Source identification 'PS3' is for PRETI supplementary list LW3 and 'PS2' for LW2. Sources with letter at 
the end, like '6a' are deblended neighbours of the source with same number in the main list. (2), (3) : source position 
for the J2000.0 equinox. (4) : source flux in microJansky in the LW2 filter with error bar or upper limit (5t w ). A 
'-' indicates that the source is not in the LW2 field. (5) : source flux in ADU/G/s in a 3" aperture. (6) : detection 
level (nr w ) in the LW2 image. (7) : source flux in microJansky in the LW3 filter with error bar or upper limit (5t w ). 
(8) : source flux in ADU/G/s in a 6" aperture. (9) : detection level (nr w ) in the LW3 image. (10) : Identificator of 
the HST field in which the source falls (hd stands for HDF, others for flanking fields,'-' mea ns ou tside of all HST 
fields). (11) : s pectro scopic redshift with source incat ed by a note, (c) : data from Cohen et al. ( |l996|) , (p) : data from 
Phillips et al. ( |l997|) , (i) : compiled in Cowie et al. (|l997|) , (h) data from Hogg et al. (|l977b|) . (12) : mo rphological 
type of the object, '(i)' indicates a type compiled in Cowie et al. ( 1997| ),'(p)' type from Phillips et al. ( |l997 ), otherwise 
determinated by us. E : Elliptical or SO galaxy. S : spiral. M : merger. P : peculiar/irregular, C : compact object from 
Phillips et al. ( 1997| ) sample, A : active nuclei galaxy, G : source falls on a group. U : unknown type. N : no optical 
counterpart. K : counterpart only in IRIM K image, st : star-like object . A '?' indicates a dubious type identification. 
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Table 5. Comparison of the LW2 catalog obtained here with PRETI, with the catalogs by IC (Goldschmidt et al. 1997) 
and IAS (Desert et al. 199? ), as well as radio sources detected at 8.4 GHz (Fomalont et al. |1997 ) 
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Table 6. Comparison of the LW3 catalog obtained here with PRETI, with the catalogs by IC (Goldschmidt et al. 1997) 
and IAS (Desert et al. 1998| ), as well as radio sources detected at 8.4 GHz (Fomalont et al. |1997 ) 
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267 


_ 


_ 


HDF_PM3_23 


144±g 


_ 


_ 


HDF_ALL_LW3_24 


130 


3646+1448 


36 


HDF_PM3.24 


307+e? 


J123648. 1+621432 


231 


HDF_ALL_LW3.25 


265 


_ 


_ 


HDF_PM3.26 


150t2 


_ 


_ 


HDF_ALL_LW3.26 


91 


_ 


_ 


HDF_PM3.27 


320+gj 


J123649. 8+621319 


471 


HDF_ALL_LW3.27 


250 


3649+1313 


22 


HDF_PM3.28 


3411*0 


_ 


_ 


HDF_ALL_LW3.28 


367 


_ 


_ 


HDF_PM3_29 


48l 3 2 9 


_ 


_ 


_ 


_ 


3651+1221 


18 


HDF_PM3.30 




J123651.5+621357S 


155 


HDF_ALL_LW3.30 


135 


_ 


_ 


HDF_PM3.31 


174lf 3 


J123653.0+621116 


327 


HDF_ALL_LW3.31 


170 


_ 


_ 


HDF_PM3.32 


180+f 3 


J123653.6+621140 


138 


HDF_ALL_LW3.32 


155 


3653+1139 


15 


HDF_PM3.33 


179+™ 






HDF_ALL_LW3.33 


143 






HDF_PS3.24 


23+i? 










3655+1311 


12 


HDF_PS3_37 


225+_™ 


J123658.1+621458S 


210 


HDF_ALL_LW3.34 


170 






HDF_PS3.39 


157+H 


J123658. 7+621212 


336 


HDF_ALL_LW3.35 


156 






HDF_PS3_40 


295l« 


J123700. 2+621455 


290 


HDF_ALL_LW3.36 


301 






HDF_PS3_32 


i5i 9 9 










3701+1147 


19 


HDF_PM3.42 


162^2 


J123702.0+621127 


326 


HDF_ALL_LW3.37 


235 






HDF_PM3.44 


1441™ 


J123702. 5+621406 


332 


HDF_ALL_LW3.38 


167 






HDF_PM3.45 


431+ 43 


J123705.7+621157 


471 


HDF_ALL_LW3.39 


359 






HDF_PS3_40 


306+^ 


J123709.8+621239 


510 


HDF_ALL_LW3.41 


267 






HDF_PS3.38 


891™ 










3708+1245 


15 



